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A  developing  surface  frontal  system  with  an  overlying  polar  front 
jet  and  subtropical  jet  stream  merging  region  was  present  over  the  con- 
tinental United  States  during  22  to  24  March,  1966.  Possible  cross- 
over flow  between  these  two  jet  streams  has  been  studied  in  detail  by 
isentropic  air  trajectories  at  four  potential  temperature  levels  from 
310  K  to  335  K, 

The  results  of  the  trajectory  analyses  on  the  four  0  surfaces 
have  indicated  a  cross— over  of  air  parcels  from  the  subtropical  jet 
stream  to  the  polar  front  jet  stream.  No  evidence  has  been  found  for  a 
reciprocal  relationship  for  air  parcels  initially  associated  with  the 
polar  front  jet  between  310  K  and  335  K  isentropic  levels. 
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SECTION  I 
INTRODUCTION 

Studies  by  Reiter  and  Whitney  (196$)  of  the  region  where  a  polar 
front  jet  and  a  subtropical  jet  stream  are  in  close  proximity  to  each 
other  have  produced  a  provocative  hypothesis  that  there  exists  a  cross- 
over flow  between  the  two  respective  jets.  One  corollary  is  that  cer- 
tain air  parcels  in  the  polar  front  jet,  rather  than  follow  a  normal 
cyclonic  trajectory  over  the  continental  United  States,  crossover  and 
proceed  along  a  trajectory  path  similar  to  that  of  air  parcels  in  the 
subtropical  jet  stream.  A  second  corollary  proposes  that,  rather  than 
follow  a  basically  anticyclonic  trajectory,  certain  air  parcels  in  the 
subtropical  jet  cross— over  and  then  continue  along  a  trajectory  path 
similar  to  that  of  air  parcels  in  the  polar  front  jet  stream  flow. 

Regardless  of  the  validity  of  the  hypothesis,  the  isentropic  tra- 
jectory of  a  representative  air  parcel  in  either  jet  can  be  computed  by 
the  method  developed  by  Danielsen  (1961).  This  computational  method 
describes  with  considerable  accuracy  the  movement  of  air  parcels  along 
an  isentropic  surface.  However,  unlike  several  other  approaches  to 
calculating  air  trajectories  which  assume  that  air  parcels  move  quasi- 
horizontally  on  a  surface  of  constant  height  or  pressure,  the  Danielsen 
approach  gives  the  necessary  third  dimension  —  a  measurement  of  the 
average  vertical  motion. 

Throughout  this  study  diabatic  effects  on  trajectory  motion  in 
this  non-precipating  region  are  neglected  because  diabatic  contribu- 
tions to  vertical  motion  are  normally  small  compared  to  adiabatic  con- 
tributions . 
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During  late  March  1966  a  strong  developing  surface  low  with  as- 
sociated frontal  systems  moved  east-northeast  across  the  continental 
United  States.  Also  associated  with  this  surface  system  were  a  strong 
polar  front  and  a  subtropical  jet  stream  in  close  proximity  to  each 
other  as  illustrated  in  Figure  1.  The  period  between  22  March  0000  GMT 
to  24  March  0000  GMT  was  selected  as  being  most  favorable  for  testing 
the  hypothesis  of  Reiter  and  Whitney  (1965). 
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SECTION  II 
NATURE  AND  ACQUISITION  OF  THE  DATA 

All  the  data  in  the  analysis  comes  from  radiosonde  and  rawinsonde 
checked  data  found  in  the  U^S.  Weather  Bureau  Northern  Hemisphere  Data 
Tabulations  for  22,  23,  and  24  March  1966. 

The  period  under  study  contains  five  observation  times,  22  0000 
GMT,  22  1200  GMT,  23  0000  GMT,  23  1200  GMT,  and  24  0000  GMT.  For  each 
observation  time,  an  area  of  interest  is  determined  by  plotting  an  area 
encompassing  not  only  the  surface  system  but  also  the  region  ot  the 
merging  jets.  Table  1  shows  the  data  areas  for  each  observation  time. 
Having  delineated  the  regions  of  interest,  the  checked  data  is  utilized 
from  each  upper-air  station  within  these  regions. 

The  data  for  each  upper-air  station  consists  of  all  standard  pres- 
sure levels  in  50  mb  intervals  from  1000  mb  to  200  mb.  From  200  mb  to 
100  mb  the  interval  is  25  mb.  Temperature  in  degrees  Celsius,  heights 
in  meters,  wind  direction  in  degrees,  and  wind  speed  in  meters  per  sec- 
ond are  tabulated  for  each  standard  level.  In  addition  to  the  standard 
levels,  the  significant  levels  are  recorded  with  their  tabulated  temp- 
eratures and  heights 
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Table  1.  Data  areas  used  for  each  observation  time. 


OBSERVATION 
TIME 

NORTH 

SOUTH 

WEST 

EAST 

22  March 
0000  GMT 

U.S.-Canada 
Border 

U.S.-Mexico 
Border  and 
Gulf  Coast 

West  Coast 
U  S 

90°W 

22  March 
1200  GMT 

U.S.-Canada 

Border 

U.S.-Mexico 
Border  and 
Gulf  Coast 

115°W 

80°W 

23  March 
0000  GMT 

U.S.-Canada 
Border 

U.S.-Mexico 
Border  and 
Gulf  Coast 

105°W 

75°W 

23  March 
1200  GMT 

U.S.-Canada 
Border 

U.S.-Mexico 
Border  and 
Gulf  Coast 

105°W 

75°W 

24  March 
0000  GMT 

L 

U.S.-Canada 
Border 

U.S.-Mexico 
Border  and 
Gulf  Coast 

105°W 

65°W 
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SECTION  III 

THE  MONTGOMERY  STREAM  FUNCTION  COMPUTER  PROGRAM 

:=  0 
Assuming  d0/dt,  the  equations  of  motion  in  6  coordinates  become 

A 


2 1;  3ixQ  2ije  3Xe 


(1) 


2v_     +CihL    +  s/£X-     =-£M     -  f u  ,  (2) 

pt  3xd         <JVd  3ye 

where  the  isentropic  stream  function  is  defined  as 

NA  =CPT+  cjZ  .  (3) 

For  the  air  trajectory  method  as  derived  by  Danielsen  (1961)  and 
later  outlined  by  Mahlman  (1965),  the  Montgomery  stream  function  is 
given  by 

K/\  =  CpTe   tgBp  f  RT(7nPP-&§)  .  (4) 

As  observed  by  Danielsen  (1959),  large  errors  of  M  result  if  val- 
ues of  Tq  and  Pq  are  selected  independently  of  each  other  from  the 
plotted  soundings „  Consequently,  for  this  program,  the  dependent  re- 
lationship of  Tq  and  Pq  on  an  isentropic  surface,  as  given  in  Poisson's 
equation,  allows  M  to  be  written 


J  r  hood  {  Tq   J 


[fOOd    V  Tq 


(5) 


as  shown  in  Mahlman  and  Kamm  (19651  M  is  now  in  a  convenient  compu- 
tational form  which  can  be  handled  readily  by  the  computer. 

The  basic  program  is  presented  in  Appendix  I  and  will  be  briefly 
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outlined  in  this  section.  After  reading  in  the  data,  all  temperatures 
are  converted  to  degrees  Kelvin.  Potential  temperatures  are  next  com- 
puted for  each  level.  For  this  study,  a  range  of  potential  temperatures 
between  290  K  and  350  K  at  5°  intervals  is  assigned  for  computation. 

The  data  s  now  searched  for  the  first  9  level  above  and  below  the 
predetermined  0  level.  With  these  two  0  values  and  the  predetermined 
0  value,  Tq  is  computed,  Pq  is  then  computed  from  Poisson's  equation 
using  Tq  and  0.  Finding  next  the  first  standard  levels  above  and  below 
the  arbitrary  isentropic  surface  and  calculating  T  by  summing  the  weight- 
ed mean  temperatures  between  all  data  points  from  the  lower  standard 
level  to  the  arbitrary  0  level,  the  Montgomery  stream  function  is  com- 
puted in  three  segments  according  to  equation  (5).  Finally,  the  wind 
direction  and  wind  speed^ifc  interpolated  according  to  the  procedure  used 
by  Mahlman  and  Kamm  (1965).  After  completing  the  calculations  for  one 
particular  isentropic  surface  at  all  stations  in  the  data  area,  the  pro- 
cedure is  repeated  for  the  next  higher  relevant  0  level. 

The  final  computer  output  gives  for  each  isentropic  surface  the  M 
value,  pressure,  static  stability,  wind  direction,  and  wind  speed  for 
all  reporting  stations.  With  these  five  parameters,  the  isentropic 
trajectories  can  be  computed  and  checked  for  accuracy,  thus  enabling 
the  three-dimensional  tracking  of  selected  air  parcels. 
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SECTION  IV 
DETERMINATION  OF  RELEVANT  ISENTROPIC  SURFACES 

Possessing  now  the  computed  parameters  M,  Pq,  wind  direction  and 
wind  speed  for  all  the  isentropic  surfaces  at  5°  intervals  between  290  K 
and  350  K,  the  next  step  in  the  procedure  is  to  determine  particular  6 
levels  which  are  representative  of  the  isentropic  surfaces  present  in 
the  polar  front  jet  and  the  subtropical  jet  streams. 

The  observation  time  23  March  12000  GMT  is  selected  as  the  synoptic 
map  most  representative  of  the  five  observation  times,  A  north-south 
cross  section  on  longitude  95°   W  is  constructed  with  plotted  values 
of  wind  speed  and  potential  temperature  as  depicted  in  Figure  2.  An- 
alysis of  this  cross  section  shows  that  four  isentropic  surfaces  might 
possibly  be  of  interest  —  310  K,  320  K,  330  K,  and  335  K.  These  four 
6  surfaces  are  chosen  because  of  their  physical  location  with  respect 
to  both  jet  streams  and  their  heights  above  sea  level,  so  that  there 
will  be  a  minimum  of  orographic  and  lower  level  interference. 

Basic  maps  are  plotted  on  1:15,000,000  scale  charts  for  the  four 
0  surfaces.  Each  of  these  isentropic  surfaces  is  analysed  for  the 
Montgomery  stream  function  using  a  0.5x10'  cm^  sec"''  interval,  for  pres- 
sure using  a  50  mb  interval,  and  for  wind  speed  using  a  10  m  sec"1  iso- 
tach  interval.  Figures  3,  4,  and  5  are  typical  examples  of  these  anal- 
yses. 


20 


fig.  *>.     23  March  1966,  1200  GMT  nortn-south  cross-section  on  longi- 
tude 95°W.  -Solid  lines  are  isotachs  (units  m  sec"1);  dashed  lines  ar*. 
potential  temperature  (°K).  310  K,  320  K,  330  K,  and  335  K  are  the 
isentropic  surfaces  chosen  for  trajectory  analysis  in  this  study. 
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SECTION  V 
COMPUTATION  OF  ISENTROPIC  AIR  TRAJECTORIES 

Mahlman's  outline  (1965)  of  the  Danielsen  method  (1961)  develops 
in  detail  the  two  working  equations  for  computation  of  isentropic  air 
trajectories, 

|^.dt  =OsA*-M,)d  f  •  (V*-V*)e  (6) 

D  =  \^Ab.  (7) 

The  left  hand  term  in  equation  (6)  represents  the  local  change  in 
the  stream  function.  The  local  change  is  graphically  determined  by  an- 
alysing the  12  hour  local  changes  of  aM  at  each  upper-air  station  or 
data  point.  For  the  four  time  periods  of  this  study,  tM   is  analysed 
using  an  increment  of  0.25x10'  cm  sec"  .  ^  j.a  then  averaged  along  the 
trajectory  path  to  obtain  an  approximation  to  the  left  hand  side  of 
equation  (6).  Figure  6  illustrates  the  local  change  of  M  for  the  period 
22  March  0000  GMT  to  22  March  1200  GMT  on  the  330  K  6  surface. 

The  second  term  on  the  right  hand  side  of  equation  (6)  expresses 
the  change  in  kinetic  energy  of  the  horizontal  wind.  After  converting 
the  units  of  this  kinetic  energy  term  to  the  units  of  M  (10?  cm  sec"*^) 
and  rei-arranging  terms,  equation  (6)  can  be  written  as 

M2=M,  +  O.00os(V,2-Vt)e+ZM  (8) 

for  V  expressed  in  m  sec""  . 
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The  units  of  equation  (7),  the  distance  equation,  are  converted  to 
degrees  latitude  for  a  6  hour  trajectory,  becoming 

DiO./?^.  (9) 

The  objective  now  is  to  locate  a  value  of  M  at  the  final  point  on 
a  trajectory  path  such  that  equations  (8)  and  (9)  are  simultaneously 
satisfied,,  Detailed  procedures  for  calculating  isentropic  air  traject- 
ories can  be  found  in  Daniels en's  paper  (1961).  However,  the  procedures 
used  for  this  study  will  be  briefly  outlined  for  one  time  period. 

For  the  first  12  hour  period,  22  0000  GMT  to  22  1200  GMT,  a  rep- 
resentative sample  of  initial  data  points  in  and  around  the  isotach 
center  of  the  polar  front  jet  and  the  subtropical  jet  streams  is  selected 
on  the  22  0000  GMT  M  analysis  for  the  four  relevant  9  surfaces.   The 
first  estimate  of  the  path  of  the  trajectory  is  obtained  by  sketching 
both  the  22  0000  GMT  and  the  22  1200  GMT  streamlines  around  one  of  the 
initial  data  points.  An  average  streamline  is  determined  by  giving 
more  weight  to  the  22  0000  GMT  streamlines  for  the  first  6  hour  period 
and  to  the  22  1200  GMT  for  the  second  6  hour  period.  A  6  hour  traject- 
ory is  computed  by  equation  (9)  using  the  wind  speed  at  the  initial  data 
point  and  is  plotted  along  the  average  streamline  using  1°  increments. 
A  mean  wind  at  the  6  hour  position  is  determined  from  the  22  0000  GMT 
and  22  1200  GMT  isotach  analyses ,  This  mean  wind  is  then  averaged 
with  the  initial  wind  to  determine  an  estimated  mean  wind  over  the 
first  6  hour  period.  Equation  (9)  is  again  employed  to  calculate  a 
better  estimate  of  the  6  hour  trajectory  position. 

Using  the  mean  wind  at  the  6  hour  position,  the  trajectory  for  the 
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next  6  hour  period  is  computed  by  equation  (9)#  The  winds  at  the  12 
hour  and  6  hour  positions  are  averaged  and  this  resulting  mean  wind  is 
used  to  re-estimate  the  12  hour  position. 

Now,  the  trajectory  computed  by  using  equation  (9)  must  be  tested 
by  equation  (8),  To,  calculate  AM,  the  trajectory  determined  by  the 
distance  equation  is  superimposed  on  the  respective  aM  analysis  and 
divided  into  four  equal  segments,  Simpson's  rule  is  utilized  to  calcu- 
late the  12  hour  local  change  of  M  integrated  along  the  trajectory  path. 
Having  aM,  V^,  V~,  and  M-^,  M^  is  computed. 

At  this  stage  in  the  procedure,  an  observed  value  of  M2,  determined 
from  the  first  estimate  to  the  trajectory  and  the  22  1200  GMT  M  analysis, 
and  a  computed  value  of  M2,  determined  by  equation  (8),  have  been  ob- 
tained, Mo-obs  ^s  then  compared  to  M2_comD.  If  their  absolute  difference 
is  equal  to  or  less  than  0.1x10'  cm2  sec  ,  the  initial  trajectory  sat- 
isfies both  equations  and  this  trajectory  represents  the  12  hour  air 
parcel  movement  on  the  isentropic  surface. 

Generally,  for  the  first  estimated  trajectory,  Mo_comr)  ""  ^2-obs 
is  greater  than  O.lxKr  cm2  sec  ,  If  this  is  the  case,  the  12  hour 
position  calculated  by  equation  (9)  is  moved  perpendicular  to  the  mean 
flow  towards  the  M2_COmp  value  obtained  from  equation  (8).  A  new  V2  at 
this  second  estimated  12  hour  position  is  determined  from  the  22  1200 
GMT  isotach  analysis. 

The  initial  trajectory  is  now  adjusted  throughout  its  entire  path 

to  conform  to  the  second  estimated  12  hour  position.  The  resulting 

differences  between  the  adjusted  trajectory  and  the  initial  trajectory 

have  changed  the  6  hour  position.  Hence,  a  new  mean  wind  over  the  first 

6  hour  period  and  the  second  6  hour  period  is  calculated  by  the  method 

outlined  previously. 

27 


Following  the  original  procedures  with  the  new  6  hour  mean  winds, 
a  second  estimated  12  hour  trajectory  position  on  the  adjusted  trajec- 
tory is  computed.  If  |M2_comp  -  M2-obs  new  is  less  than  or  e^11*1  to 
0.1x10'  cm^  sec~2,  the  second  estimated  trajectory  satisfies  the  dis- 
tance and  the  energy  equations  simultaneously.  If  not,  the  previous 
procedures  are  again  repeated  until  convergence  is  reached.  Thus,  an 
air  parcel  has  been  successfully  tracked  for  12  hours  on  an  isentropic 
surface  by  successive  approximations. 

Figure  7  illustrates  the  final  estimated  trajectories  for  the 
22  0000  GMT  to  23  0000  GMT  time  period  on  the  330  K  isentropic  surface. 
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SECTION  VI 
RESULTS  OF  THE  ISENTROPIC  AIR  TRAJECTORY  COMPUTATIONS 

After  computing  a  representative  sample  of  trajectories  for  each 
of  the  four  isentropic  surfaces,  (310  K,  320  K,  330  K,  and  335  K),  for 
all  four  time  periods,  the  trajectories  are  superimposed  on  the  respec- 
tive pressure  analyses,  so  that  Pq  can  be  determined  for  the  initial  and 
final  position.  Ctf  is  then  calculated  for  each  12  hour  trajectory,  where 

To  illustrate  the  three-dimensional  air  parcel  movement,  an  analy- 
sis of  trajectory  position  relative  to  the  surface  front  position  is 
made.  The  positional  relationship  between  a  trajectory  and  the  moving 
surface  frontal  system,  as  shown  in  Figure  7,  appears  quite  consistent 
for  all  four  0  surfaces.   Trajectory  2  on  the  330  K  surface,  initially 
over  Las  Vegas,  Nevada  (Station  336),  has  an  average  downward  vertical 
motion  (65  positive)  of  +16  mb/12hrs  for  the  22  0000  GMT  to  22  1200  GMT 
period.  This  downward  vertical  motion  is  reasonable  since  the  Trajec- 
tory 2  position  is  generally  lagging  the  cold  front  position  for  the 
12  hour  period  —  thus  remaining  in  a  subsidence  region. 

As  Trajectory  2  progresses  to  the  northeast,  its  position  relative 
to  the  frontal  system  is,  on  the  average,  leading  the  surface  system. 
This  leading  situation  is  indicative  of  the  average  upward  vertical 
motion  of  -1  mb/12hrs  for  the  22  1200  GMT  to  23  0000  GMT  period. 

This  leading  and  lagging  relative  relationship  between  isentropic 
air  trajectories  and  the  surface  frontal  system  was  consistent  for  all 
time  periods  and  isentropic  surfaces.  Thus,  through  the  laborious  4 
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level  computations,  a  physical  understanding  of  air  trajectory  three- 
dimensional  motion  has  been  attained.  In  addition,  a  mental  picture  of 
this  motion  can  be  generated  by  relating  these  trajectories  to  the  sur- 
face phenomenon. 
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SECTION  VII 
TESTING  THE  HYPOTHESIS 

Recent  studies  by  Danielsen  (1961),  Reiter  and  Whitney  (1965),  and 
Mahlman  (1965)  have  investigated  cases  of  apparent  jet  stream  cross-over 
at  isentropic  levels  below  310  K.  Assuming  the  cross-over  phenomenon 
does  exist  at  the  lower  levels  of  the  atmosphere,  it  is  reasonable  to 
inquire  as  to  the  existence  of  a  compensatory  crossover  flow  above 
the  310  K  isentropic  level.  Therefore,  this  hypothesis  will  be  tested 
by  further  analyses  of  trajectories  associated  with  the  respective  jet 
streams  within  the  merging  region. 

The  isentropic  air  trajectories  are  analysed  for  their  positions 
relative  to  isotach  centers  of  the  polar  front  jet  and  the  subtropical 
jet  streams.  The  trajectories  are  further  analysed  for  their  positions 
relative  to  the  region  of  apparent  jet  stream  merging.  As  illustrated 
in  Figure  3,  Trajectory  1  at  22  0000  GMT  is  located  within  the  50  m 
sec"  isotach  center  of  the  polar  front  jet  stream  on  the  335  K  sur- 
face. Maintaining  a  cyclonic  curvature  with  an  average  downward  vertical 
motion  of  +8  mb/12hrs  for  the  succeeding  12  hour  period,  an  air  parcel 
associated  with  Trajectory  1  remains  in  the  polar  jet  flow  and  by  22 
1200  GMT  is  well  out  of  the  merging  region.  Referring  again  to  Figure  8, 
an  air  parcel  associated  with  Trajectory  2  follows  a  similar  path  as 
that  of  Trajectory  1,  but  is  displaced  more  to  the  south  throughout 
the  22  0000  GMT  to  22  1200  GMT  period.  As  before,  this  trajectory  can 
be  tracked  through  the  merging  region  and  on  to  thj  northeast. 

After  analysing  the  trajectories  associated  with  the  polar  front 
jet  stream  on  the  four  9  surfaces  for  the  four  time  periods,  the  mean 
three-dimensional  flow  of  air  parcels  in  the  polar  front  jet  is  deter- 
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mined  to  be  basically  cyclonic  in  nature  —  descending  when  moving  south- 
east and  ascending  when  moving  northeast.  This  mean  flow  corroborates 
previous  studies. 

For  the  22  0000  GMT  to  22  1200  GMT  period,  the  apparent  jet  merging 
region  is  over  New  Mexico  and  Arizona.  Trajectory  3  of  Figure  8  at 
22  0000  GMT  begins  within  the  merging  region  over  Tucson,  Arizona 
(Station  274).  An  air  parcel  following  this  trajectory  appears  to 
disassociate  from  the  main  flow  of  the  subtropical  jet  and  proceeds  to 
the  northeast  with  an  average  upward  vertical  motion  of  -10  mb/12hrs 
for  the  first  12  hour  period.  By  22  1200  GMT,  the  air  parcel  associated 
with  Trajectory  3  has  entered  the  isotach  center  over  Missouri  and 
continues  to  move  northeastward  toward  the  Great  Lakes  region.  This 
phenomenon  of  apparent  disassociation  of  air  parcels  from  the  subtropical 
stream  flow  is  evident  on  the  four  0  surfaces  for  the  trajectories  with 
initial  position  over  Tuscon  at  22  0000  GMT. 

Trajectory  4  of  Figure  8  is  located  over  El  Paso,  Texas  (Station 
270)  at  22  0000  GMT0  Unlike  Trajectory  3,  an  air  parcel  associated  with 
Trajectory  4  remains  within  the  anticyclonically  curving  subtropical 
jet  stream  flow  for  the  succeeding  24  hour  time  period. 

To  substantiate  the  statement  made  in  the  previous  paragraph  con- 
cerning apparent  disassociation  of  air  parcels  from  the  subtropical  jet 
stream  flow  on  a  9  surface  other  than  335  K,  refer  again  to  Figure  6. 
Trajectory  3  on  the  330  K  surface  also  begins  over  Tucson,  Arizona  at 
22  0000  GMT0  This  trajectory  path  appears  to  cross-over  and  by  22  1200 
GMT  possesses  the  characteristic  flow  pattern  of  the  polar  front  jet 
stream. 

By  23  March  0000  GMT,  the  jet  stream  merging  region  is  located  over 

Texas,  Trajectory  2  on  the  330  K  isentropic  surface  in  Figure  9, 
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initially  located  over  Del  Rio,  Texas  (Station  261),  appears  to  continue 
within  the  subtropical  jet  influence  and  maintains  an  average  upward 
vertical  motion  of  -29  mb/12hrs  for  the  23  0000  GMT  to  23  1200  GMT 
time  period. 

The  isotach  analyses  of  the  four  isentropic  surfaces  for  the  23 
1200  GMT  to  24  0000  GMT  period  define*  the  merging  region  of  the  polar 
front  and  the  subtropical  jet  streams  most  clearly.  Trajectory  3  of 
Figure  10  was  selected  because  not  only  was  the  air  parcel  following 
this  trajectory  path  associated  with  the  subtropical  jet  flow,  but  the 
23   1200  GMT  position  over  Midland,  Texas  (Station  265)  was  well  within 
the  merging  region.  Trajectory  3  on  the  320  K  and  330  K  isentropic 
surfaces  appears  to  crossover  and  proceed  cyclonically  to  the  north- 
east, apparently  identifying  itself  with  the  polar  front  jet  flow  at 
these  levels  as  illustrated  in  Figure  10.  However,  trajectories  on  the 
335  K  and  310  K  surfaces  located  over  Midland  at  23  1200  GMT  are  dis- 
placed to  the  south  of  the  320  K  and  330  K  trajectory  position,  indicating 
a  continued  association  with  the  subtropical  jet  stream  flow.  It  is 
interesting  to  note  that  Trajectory  k,   located  initially  on  the  anti- 
cyclonic  side  of  the  subtropical  jet  stream,  apparently  remains  well 
within  this  flow. 

The  subtropical  jet  is  located  north  of  its  normal  position  for  the 
23  March  1200  GMT  to  24  March  0000  GMT  period.  Trajectory  3  of  Figure 
10  is  initially  located  on  the  cyclonic  side  of  the  jet  and  in  the 
central  section  of  the  merging  region.  Trajectories  computed  to  the 
south,  or  on  the  anticyclonic  side  of  the  subtropical  jet  stream,  follow 
the  classical  anticyclonic  path  over  the  southern  United  States  as 
shown  in  Figures  7  and  9. 
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Potential  vorticity  (P)  of  an  air  parcel  is  known  to  be  a  conserva- 
tive property  for  adiabatic  and  frictionless  air  motions.  As  a  further 
test  for  the  four  trajectories  with  apparent  cross-rover  flow,  P  is  com- 
puted for  the  initial,  12  hour,  and  24  hour  isentropic  trajectory  posi- 
tions. If  the  three  P  values  do  not  vary  appreciably,  this  lends  further 
credence  to  the  validity  of  the  trajectories  indicating  crossfover  flow. 
Consequently,  the  fact  that  potential  vorticity  is  conserved,  as  shown 
in  Table  2,  helps  substantiate  the  results  of  the  isentropic  trajectory 
computations. 

As  shown  in  Figure  8,  Trajectory  3  on  the  335  K  isentropic  surface 
for  the  first  12  hours  follows  a  cyclonic  path  and  for  the  second  12 
hours  an  anticyclonic  path.  If  the  assumption  of  conservation  of  po- 
tential vorticity  is  valid  for  this  assumed  frictionless  and  adiabatic 
motion,  then  the  large  24  hour  relative  vorticity  change  which  reflects 
this  curvature  change  must  be  explained  by  a  significant  increase  of 
one  of  the  parameters  in  the  P  equation.  One  possibility  is  a  rapid 
increase  in  the  Coriolis  parameter,  but  for  this  trajectory,  variations 
in  f  are  small  as  shown  in  Table  2,  However,  Table  2  shows  that  the 
relative  vorticity  decrease  is  compensated  by  a  rapid  static  stability 
increase.  To  some  extent  this  argument  helps  justify  the  physical 
mechanisms  involved  in  the  crossover  phenomenon. 
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Table  2.  Approximate  values  of  potential  vorticity  computed  at  rela- 
vant  trajectory  positions.  Note  the  static  stability  and  absolute  vor- 
ticity values  corresponding  to  respective  trajectory  curvatures  shown 
in  Figures  7,  8,  and  10. 


9  SURFACE 


330  K 

335  K 

320  K           !       330  K 

FIGURE 
TRAJECTORY 

7 
3 

8 
3 

10 
3 

10 
3 

DATE 

22  0000  GMT 

22  1200  GMT 

23  0000  GMT 

22  0000  GMT 

22  1200  GMT 

23  0000  GMT 

23  1200  GMT 

24  0000  GMT 

23  1200  GMT 

24  0000  GMT 

STATIC 
STABILITY 

(deg  cm  sec     g     ) 

3. 8x10"  5 
4. 5xl0"5 
8. 0x10* 5 

3.8xl0"5 
5.4x10*5 

8.0x10*5 

0.53xlO~5 
0.50x10*5 

5.6x10*5 
5.3x10*5 

CORIOLIS 

PARAMETER 

(sec"~l) 

0.73xl0~4 
0.83xl0~4 
0.93xlO~4 

0.73x10*4 
0,83xl0*4 
0.93x10*4 

0.79x10*4 
0.87x10*4 

0.79x10*4 
0.83x10*4 

SHEARING 
VORTICITY 

(sec"-1-) 

0.30x10*4 
0.0 
0.0 

0.30x10*4 

0.12x10*4 

-0.06xl0~4 

-0.20xl0"4 
-0.25x10*4 

0.0 
-O.OoxlO"4 

CURVATURE 
VORTICITY 

(sec-1) 

0.28x10*4 

0.27x10*4 

-0.28x10*4 

0.27x10*4 

0.31x10*4 

-0.22x10"^ 

0.27x10*4 
0.34xl0~4 

0.32x10*4 
0.30xl0"4 

ABSOLUTE 
VORTICITY 

(sec"1) 

1.31x10*4 
l.lOxlO*4 
0o65xl0~4 

1.30x10*4 
1.10x10*4 
0.65xlO*4 

- 
0.86x10*4 

0.96x10*4 

l.llxlO"4 
1.19x10*4 

POTENTIAL 
VORTICITY 

(deg  cm  sec     g*1) 

5xl0~9 
5xl0~9 
5x10*9 

5xlO*"9 
6x10-9 
5x10*9 

5x10*9 
5x10*9 

6x10*9 
6x10-9 
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SECTION  VIII 
CONCLUSIONS 

As  indicated  in  Figures  7 ,   8,  and  10  and  their  subsequent  analyses, 
cross-over  of  air  parcels  in  the  subtropical  jet  stream  flow  to  the 
polar  front  jet  flow  apparently  exists  between  the  310  K  and  335  K 
isentropic  surfaces.  This  conclusion  is  dependent  upon  the  validity  of 
the  adiabatic  assumption.  However, refer ring  again  to  Figures  7,8 
and  10,  it  can  be  observed  that  small  directional  wind  shears  in  the 
vertical  exist  through  the  10°  isentropic  layer.   Therefore,  diabatic 
effects,  even  if  present,  would  not  appreciably  affect  the  results.  An 
additional  implication  from  these  analyses  is  the  existence  of  a  mass 
transport  phenomencn  in  the  jet  stream  region  as  observed  by  Reiter  and 
Whitney  (1965). 

The  cross-over  phenomenon  of  air  parcels  initially  in  the  polar 
front  jet  stream  is  not  observed  in  this  case  study  for  the  relevant 
isentropic  surfaces.  However,  the  apparent  lack  of  air  parcel  cross- 
over bo  the  subtropical  jet  stream  flow  at  these  0  levels  might  be  ex- 
plained by  the  limited  number  of  isentropic  surfaces,  synoptic  situa- 
tions, and  parameters  utilized  in  this  research. 
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APPENDIX  I 


The  Montgomery  stream  function  computer  program  developed  for  this 
study  (Sections  A  and  B)# 
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APPENDIX  I;  Section  A,  List  of  symbols  used  in  the  computer  program. 


SYMBOL 

DD1 

DD2 

DDTH 

F 

Fl 

F2 

F3 

FF1 

FF2 

FFTH 

PP 

PP2 

PPB 

PPT 

PTH 

STAB 

T 

TB 

TBAR 

TH 

THB 

TH1 

THT 

TT 


MEANING 
wind  direction  at  first  standard  level  below  TH1 
wind  direction  at  first  standard  level  above  TH1 
interpolated  wind  direction  at  TH1 
Montgomery  stream  function 
first  segment  of  F 
second  segment  of  F 
third  segment  of  F 

wind  speed  at  first  standard  level  below  TH1 
wind  speed  at  first  standard  level  above  TH1 
interpolated  wind  speed  at  TH1 
pressure  at  first  standard  level  below  TH1 
pressure  at  first  standard  level  above  TH1 
pressure  at  first  level  below  TH1 
pressure  at  first  level  above  TH1 
pressure  at  TH1 


stability  parameter  =  - 


a©. 

9P 


temperature 

T  at  first  level  below  TH1 

average  T  of  layer  between  PP  and  PP2 

potential  temperature 

TH  at  first  level  below  TH1 

assigned  TH  value 

TH  at  first  level'  above  TH1 

T  at  first  level  above  TH1 
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SIMBOL  MEANING 

TTH  T  at  THl 

ZTH  geopotential  height  of  TH1 
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APPENDIX  I:     Section  B#     The  Montgomery  stream  function  computer  program. 


MONTGOMERY    STREAM    FUNCTION    COMPUTER    PROGRAM    —    I    JUNE    1968 

RADIOSONDE    DATA    FOR      22    MARCH    0000    GMT    TO    24   MARCH    0000    GMT 
1966.       FIVE    OBSERVATION    TIMES    —    RESPECTIVE    MATRIX    A,B,C,D,E. 

DIMENSION  A(1190,5),  B(1478,5),  0(1290,5),  E(1484,5),  C(1270,5) 
1  G(1484,5),  TH(1484).  T(1484) 

READ  (5, 5001  ((A(I,J),  J  =  1,5).  I  =  1,  1190) 

READ  (5,500)  ((B(I,J),  J  =  1,  5),  I  =  I,  1478) 

READ  (5,500)  ((C(I,J),  J  »  1,5),  I  =  1,  1270) 

READ  (5,500)  ((D(I,J),  J  *  1,5),  I  =  1,  1290) 

READ  (5.500)  ((Ed, J),  J  *  1,5),  I  =  1,1484) 
500  FORMAT  (10F8.0) 

WRITE  (6,600) 
600  FORMAT  (IH1) 

N  =  1 

1  IF  (N  .EO.  1)  GO  TO  2 
IF  (N  .EQ.  2)  GO  TO  3 
IF  (N  .EQ.  3)  GO  TO  4 
IF  (N  .EQ.  4)  GO  TO  5 
IF  (N  .EQ.  5)  GO  TO  6 
IF  (N  .EQ.  6)  GO  TO  17 

2  M  =  1190 
GO  TO  7 

3  M  =  1478 
GO  TO  7 

4  M  =   1270 
GO  TO  7 

5  M  =  1290 
GO  TO  7 

6  M  =  1484 

7  00  10  I  =  It  M 
DO  10  J  =  1.  5 

IF  (N  .EQ.  1)  GO  TO  20 
IF  (N  .EQ.  21  GO  TO  21 
IF  (N  .EQ.  3)  GO  TO  22 
IF  (N  .EQ.  4)  GO  TO  23 
IF  (N  .EQ.  5)  GO  TO  24 

20  G(I,J)  =  A( I, J) 
GO  TO  10 

21  G(I,J)  -  B( I, J) 
GO  TO  10 

22  G(ItJ)  =  C( I, J) 
GO  TO  10 
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23  G( I.J)  -  D( I, J) 
GO  TO  10 

24  G( I.J)  ■  E( I, J) 

10  CONTINUE 

25  00  11  I  *  1.  N 

11  WRITE  (6.610)  I.  (G(I,J),  J  =  1,  5) 
610  FORMAT  (20X,  110,  5X,  5F10.1) 

C 

C      CONVERT  TEMPERATURES  AND  COMPUTE  POTENTIAL  TEMPERATURES 


C 


c 


WRITE  (6,612) 
612  FORMAT  ( 1H1) 


DO  12  I  =  1,  M 
IF  (GU.3)  .EO.  777.0)  GO  TO  8 
T( I)  =  6(1,2)  ♦  273.2 

TH(I)  =  T(I)  *  ((1000.0  /  G(I,1»)  **  0.2857) 
GO  TO  18 
e  T( I)  =  777.0 
TH(I )  =  777.0 
18  WRITE  (6,611)  I,  T(I),  TH(I) 
611  FORMAT  (5X,  110,  2F10.1) 
12  CONTINUE 
C 
C      CHECK  FOR  DISCREPENCIES  IN  POTENTIAL  TEMPERATURES  ABOVE  800  MB 


WRITE  (6,620) 
620  FORMAT  (1H1,  T10,  'DISCREPENCIES  IN  THETA  VALUES1,  ///) 
DO  13  I  *  1.  M 

IF  (TH(I)  .EO.  777.0)  GO  TO  13 
IF  (G(I,  1)  -  800.0)  40,  40,  13 

40  IF(THd-H)  -  TH(D)  41,  13.  13 

41  WRITE  (6.630)  I,  TH(I).  THU  +  1) 
63C  FORMAT  (  5X ,  HO,  2F10.1,  ///  ) 

13  CONTINUE 

ASSIGN  POTENTIAL  TEMPERATURES  FOR  COMPUTATION 

TH1  =  290.0 
GO  TO  50 

14  TH1  =  TH1  ♦  5.0 

IF  (TH1  .GT.  350.0)  GO  TO  16 
50  WRITE  (6,  640) 
640  FORMAT  ( 1H1) 
FFTH  =  0.0 
DDTH  =  0.0 
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SEARCH  DATA  FOR  VALUES  ABOVE  AND  BELOW  LEVEL  BEING  COMPUTED 

DO  15  I  *  1,  M 

IF  (TH(I)  .EO.  777.0)  GO  TO  15 
IF  (TH(I»  -  TH1)  51,  52,  53 

51  GO  TO  15 

52  TTH  =  T(  I) 
GO  TO  54 

53  THT  ■  TH(I ) 
TT  =  T(I  ) 
PPT  =  G(Iil) 

IF  (TH(I-l)  .EO.  777.0)  GO  TO  15 
IF  (TH(I-l)  -  TH1)  55,  15,  15 
55  TB  =  T(I-l) 

PPB  =  G(I-l.l) 
THB  =  TH(I-l) 

COMPUTE  PRESSURE  AND  TEMPERATURE  AT  LEVEL  OF  INTEREST 

TTH  =  TB  ♦  (TT  -  TB)  *  ((TH1  -  THB)  /  (THT  -  THB)) 

54  PTH  =  1000.0  *  ((TTH  /  TH1 )  **  3.5001) 

FIND  NEAREST  STANDARD  LEVEL 

IPP  =  PTH  /  50.0 
AAA  =  IPP  *  50 
L  =  I 

60  IF  (G(L,1)  .EO.  AAA)  GO  TO  61 
L  =  L  ♦  1 

GO  TO  60 

61  PP  =  G(L,1)  +    50.0 
J  =  L 

62  IF  (G(J,1)  .EO.  PP)  GO  TO  63 
J  =  J  -  1 

GO  TO  62 

63  PP  =  G(J,1) 
HGT  =  G(J,3) 
FF1  =  G(J,5) 
DD1  =  G(J,4) 
AA  =  0.0 

70  IF  (PTH  -  G(J+1,D)  71,  72,  72 

71  AA  *  AA  +  (T(J)  ♦  T(J  +  1)>  *  (G(J,1)  -  G(J  +  1,1M 
J  =  J  ♦  1 

GO  TO  70 
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72  AA  =  AA  ♦  (T(J)  ♦  TTH)  *  (G(J,1)  -  PTH) 
74  TBAR  =  (1.0  /  12. G  *  ( PP  -  PTH)))  *  AA 

PP2  =  PP  -  50.0 
6*  IF  (G(J,i)  ,E0.  9?2)    GO  TO  65 

J  =  J  ♦  1 

GO  TO  64 
65  PP2  =  G( J,l) 

FF2  =  G( J, 5) 

002  =  G(J,4) 
C 

C      COMPUTE  MONTGOMERY  STREAM  FUNCTION  IN  THREE  SEGMENTS 
C 

Z  =  1.E6 

Fi  =  10.046  *  TTH  *  Z 

F2  =  9.8  *  HGT  *  1.E4 

F3  =  2.8704  *  I    *  T8AR  *  (ALOG(PP)  -  ALOG(PTH)) 
C 

C      INTERPOLATE  WIND  OIRECTION  AND  WIND  SPEED 
C 

IF     (FF1     .LT.     888.0)    GO    TO    80 

FFTH    =    888.0 

DOTH  =  888.0 
8C  IF  (FF2  .LT.  888.0)  GO  TO  81 

FFTH  =  888.0 

DOTH  =  888.0 

GO  TO  88 

81  FFTH  =  FF1  ♦  (FF2  -  FF1)  *  [[PP    -  PTH)  /  50.0) 
D03  =  DD2  -  DD1 

IF  (DD3  -  180.0)  82,  83,  83 
8  3  DO 3  =  DO 3  -  360. 0 
GO  TO  85 

82  IF  (0D3*  180.0)  84,  85,  85 

84  003  =  DD3  ♦  360.0 

85  DOTH  =  DD1  ♦  DD3  *  ( < PP    -  PTH)  /  50.0) 
IF  (DDTH)  86,  88,  87 

86  DOTH  =  DDTH  ♦  360.0 
GO  TO  88 

P7  IF  (DDTH  -  360.0)  83,  89,  89 

89  DOTH  =  DDTH  -  360.0 

88  CONTINUE 
C 
C      COMPUTE  STABILITY  FACTOR 


C 


STAB  =  (THT  -  tho)  /  ( PPB  -  PPT) 
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C      COMPUTE  MONTGOMERY  STREAM  FUNCTION 
C 

F  =  Fl  ♦  F2  «•  F3 
C 
C      COMPUTE  HEIGHT  OF  THETA  SURFACE 


C 


ZTH  =  HGT  -  (2.8704  *  TBAR  /  9.8)  *  (ALOG(PTH)  -  ALOG(PP)) 

WRITE  (6,650)  I,  STAB,  Fl,  F2,  F3,  F 
650  FORMAT  (5X,  HO,  5X,  F6.3,  5X,  3P4E15.5) 

WRITE  (6,660)  TH1,  PTH,  TTH,  ZTH,  DOTH,  FFTH 
660  FORMAT  (5X,  6F10.lt  ///) 

15  CONTINUE 
GO  TO  14 

16  N  =  N  ♦  1 
GO  TO   1 

17  STOP 
ENO 
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